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One dimensional systems sometimes show pathologically slow decay of currents. This robustness 
can be traced to the fact that an integrable model is nearby in parameter space. In integrable models 
some part of the current can be conserved, explaining this slow decay. Unfortunately, although this 
conservation law is formally anticipated, in practice it has been difficult to find in concrete cases, 
such as the Heisenberg model. We investigate this issue both analytically and numerically and 
find that the appropriate conservation law can be a non-analytic combination of the known local 
conservation laws and hence is invisible to elementary assumptions. 
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I. INTRODUCTION 

Although integrable systems are well understood in 
classical physics, their quantum counterparts are much 
more mysterious. In classical physics we have phase- 
space, and an integrable system with 2A^-dimensional 
phase-space can be characterised by a decomposition into 
a set of independent A^-dimensional tori on which the 
trajectories are constrained, circling each 'handle' inde- 
pendently with a characteristic frequency This re- 
duction in the motion of the trajectories is accomplished 
by conservation laws, of which there are N, the values 
of which can be used to annotate the distinct tori. In 
quantum physics no such generic picture exists, although 
some commonality between systems is emerging. The 
majority of integrable systems are solved using the Bethe 
Ansatz[2] which has provided a collection of 'non-trivial' 
systems which have been widely studied Q . As well as the 
spectral properties, in a parallel investigation the conser- 
vation laws have been found: Indeed, a countable num- 
ber of conservation laws can be generated by a 'boost or 
ladder operator' where the subsequent conservation law 
can be generated via commutation of the boost operator 
with the previous lawQ. Integrable systems have quite 
surprising properties. 

Various (more or less obscure but quite interesting) 
anomalies have been discovered. The additional conser- 
vation laws mean that nearby eigenstates behave inde- 
pendently and the local eigenvalue statistics is seen to 
be Poisson and not Wigner-DysonQ as it is in normal 
systems. Indeed, this property can be used as a test for 
integrabilitylQ If a current is started in an integrable 
system then it does not decay away completely, even if 
you wait forever and even at finite temperature! This 
final 'anomaly' is even of experimental interest, because 
anomalously slow decay of currents has been observed in 
systems which have an integrable model nearby in pa- 
rameter space. 



It is this final issue of long-time residuals that we ex- 
amine here, asking the question of whether the long-time 
residual is saturated by the standard conservation lawsQ 
or whether there might be a secondary mechanism or 
perhaps even 'non-local conservation laws' Q which con- 
tribute. Some aspects of these issues have previously 
been addressed but have led to the current confusion: 
It has been shown that formally the long-time residuals 
are necessarily interpretable as being due to conservation 
lawsQ and yet for a par ticular system with a 'known' 
non-vanishing residual|lOl. ITl| all the conservation laws 
yield a vanishing contribution to this residual due to 
symmetry [12[. This conundrum is our target. We cast 



doubt on the previous analysis by showing that it is not 
enough to take a complete set of mutually commuting 
symmetries to describe the residual, but then find that 
this problem is not encountered for the concrete example 
that we study in depth. 

Our investigation comes in two pieces. Firstly we ex- 
amine the problem formally, re-deriving the original sat- 
uration idea and extending it to the presence of 'mutually 
commuting' conservation laws and to the issue of compat- 
ibility of those laws with the current. Secondly we exam- 
ine a particular case numerically, showing how our ideas 
work in practice and shedding light on the previously 
mentioned conundrum. We find it possible to saturate 
the residual, but the conservation law required is non- 
analytic in the natural local conservation laws prevalent 
in the literature. 



II. FORMAL ANALYSIS 

The physical idea behind a long-time residual is that 
at some time, say t=0, a current is started and then at 
some much later time the question of how much of the 
current is still flowing is addressed. If a finite fraction 
of the current remains flowing for all time then there 
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is a long-time residual. This issue may be investigated 
through a calculation of 

(j(oo)j(O)) ^ lim 1 r rft(j(t)i-(0)), (1) 
where we employ the Heisenberg representation 

d{t) = e'"'de-'"' (2) 
and the angled brackets denote a thermal average 
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(3) 



Note that there is a linear response aspect to this issue 
and that an infinitesimal current is actually introduced 
and then propagates. The time-average is a technical 
trick to facilitate the calculation, allowing fluctuating 
terms to be eradicated: Only terms with an essentially 
static phase can contribute, others averaging to zero. 

Initially we will work with a finite system at the for- 
mal level in order to clarify the relationship between this 
quantity and conservation laws. Wc start out by using an 
eigenbasis which diagonalises the Hamiltonian, denoted 
by I n, m) and satisfying 

H \n,m) ^ e„ | n,m). (4) 



where we employed the idea that the eigenspectrum was 
discrete and we have not taken the thermodynamic limit. 



and, if wc assume that j is Hermitian, the matrix ele- 
ments between states of differing energy are lost in the 
long-time limit. To achieve a long-time residual the cur- 
rent operator must have a macroscopic component diag- 
onal in energy. 

The next important issue is that of conservation laws 



The label n denotes the distinct eigenvalues and the label 
TO is currently an arbitrary degeneracy label. As we shall 
see, degeneracy is the crucial element to the problem and 
how we deal with this constitutes a formal resolution of 
the issue. In this basis we find that 

n 771 



with the partition function 

^-E^^'^Ei-E^^-e""^"' (6) 



where we have introduced the degeneracy of the energy 
e„, as Mm and the matrix element 

0™y^(n,TO|6|n',m'). (7) 

The long-time residual may now be calculated using this 
basis and we find 



(8) 

I 

This quantity is to be compared with the initial current 



(9) 

I 

and at its most basic this constitutes operators which 
commute with the Hamiltonian. Employing the eigenba- 
sis to represent the general operator 

^ = EE (10) 

Tin' mm' 



(jMj(O)) ^^linr^^ rrftl^e"^'"5^(7.,m | e-"*^ J] e'^"'* | n', m')(n', to' | j | n,m) 

^ n m n' m' 
\^ \^ -f3en -mm' -rn'm i- _L J_ _ \^ p-/3e,i -mm' -m m 



7in' mm' 



(j(O)j(O)) = l^e-^^" ^(n,TO I J I r.',TO')(n',m' | j \ n,m) ^ | ^ E ^'w 'j,?!"^, 
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and recognising the Hamiltonian as 

H = ^^\ n,m)e,i{n,m |, 



(11) 



we find that 

^1 = E E I ^' - ^"'] Onrf in', m' I, (12) 



nn' 7nm' 



and so any operator which involves a unique energy pro- 
vides a conservation law. Indeed, any operator provides 
a conserved piece 



(13) 



and we can decompose any operator into a conserved and 
non-conserved piece 



= 0c + Oc- 



(14) 



Employin g t his decomposition we can verify Mazur's 
inequality [l3| and the saturation theorem Q 

(i(«))S(0)) = (ie(^)4(0)) + (ie(oo)4(0)) 

+ (ig(oo)B,(0)) + (ie(^)Se(O)). (15) 

Now since the conserved operators commute with the 
Hamiltonian, Oc{t) = Oc(0) = Oc, wc find that 



and the first term reproduces the previous long-time 
residual expression when A = j = B. The second and 
third terms vanish because they are off-diagonal and the 
final term vanishes because it involves non-trivial time 
dependence. We find that 

(j(oo)j(O)) = ^^^")y^f("^ (17) 

{jcjc} 

and further, for a set of orthogonal conservation laws Ca 
which satisfy 



(18) 



then 



0(oo)j(o))>E '^'°'f:!!f!'''°» . (M) 

which is Mazur's inequality. This can be proved by rep- 
resenting 



E "aOa -I- Ajc 



with 



{cfjc) {cfjm 



(CaCa) {CaCa) 



(20) 



(21) 



= ^ E E ^nn'B:::^, (16) and then 



n mm' 



(j(^)j(o))-E 



(j(0)C'a)((7tj(0)) 

(ClCa) 



{jcJc 
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The identity is the dominant content of the argument 
by Suzuki [0] and indeed shows that long-time residues 
may be interpreted as stemming from conservation laws 
in this broad sense. 

In intcgrable systems there are a collection of non- 
trivial conservation laws that are thought to play a simi- 
lar role to those in classical systems. The current issue is 
whether Mazur's inequality [l3| can be exhaustedhy these 
laws and the conserved quantity, jc, can be written down 
as a function of these conservation laws. 

The initial issue is as to which conserved operators 
are accessible to these conservation laws. We expect to 
be dealing with mutually commuting operators which in- 
clude the Hamiltonian 



0, 



(23) 



with Co — H . Since any collection of commuting Her- 



mitian operators may be simultaneously diagonalised, we 
are led to a restricted basis, | n, c, to), which satisfies 



Ca I n,c,m) = e^^ | n,c,m), 



(24) 



where for each choice of the distinct pair {n, c} and 
{n', c'} we have an a for which e^^ ^ e°,^, and the states 
are non-degenerate in this generalised sense. Now clearly 

= ^ I n,c,m)e'^^{n,c,m \=Y,<cPnc, (25) 



and each operator can be constructed from the projection 
operators 



Pnc = \ n,c,m) {n, c, m 



(26) 
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Further, it should be clear that 



n n 

n' c' ' " 



(27) 



for all of the operators described by set A. So, each 
projection operator can be created from the conservation 
laws and each conservation law can be created from the 
projection operators and the two operator subspaces are 
equivalent. The most general conservation law generated 
by this class of conservation laws takes the form 



Ac = ^ ^ I n, c, ■m)Anc{n, c, m 



A„cP«c. (28) 



nc m. 



Although the long-time residual may be generated 
from the conservation law jc, any particular set of con- 
servation laws may only generate part of the residual. 
For the subspace annotated by nc, we can only find the 
piece of jc which is proportional to the identity operator 
and consequently the maximal contribution to the Mazur 
inequality fEg fTO)) from a set of conservation laws comes 
from the conservation law 



JM 



1 



^jr;„™'(n,c,m|, (29) 



where Mnc is the degeneracy of the appropriate subspace 
and the maximal Mazur contribution is 



(iMi(0))..4^e-^-^^, 



m m -m m 
nc ncJnc nc 



(30) 



If we only employ the Hamiltonian itself as the system of 
conservation laws, then we also find 



(j (oo)3-(o))^ ^^Y. E E J- Tcjr; : 

n, rr' 



c' I 



cc' mm' 



whereas the full long-time residual in this basis is 



0-(c»)j(0))^^^ 



\ ^ p-/3en \ ^ „-m m' -m' m 

/ J Jncnc'Jnc' nc 



We find a hierarchy of four quantities: 

/o = (mm 

^ / J ^ / ^ / J Jnc n'd Jn'd ' 

n cm n'c'm' 

ry / J ^ / J / J Jnc nc' Jnc' net 



cm c'm' 



m m -m m 
nc ncJnc nc ' 



(31) 



(32) 



(33) 



(34) 



(35) 



^3 = 6-(^)j(0))ff 

_ i_ \ ^ p-/3cre \ " ^ \ ^ rj^n m-m' m' /oc\ 

^ / y "-• / y / J Jnc ncJnc' nc' ■> \'^'^ J 

n cm ^ c'm' 

with Jo > Ii > I2 > Is- Note that for a Hermitian 
matrix, the identity 



( Ha'a - 5aa' ^ Hhb j ( H aa' — ^aa' ^ ^ Hbb j — ^ Ha'aHaa' — ^ ^ HaaHa'a' > 
aa' \ b / \ b / aa' aa' 



(37) 



r 



completes the proof of the inequalities. 

Iq is the original current and the components of cur- 
rent which connect states with different energies are lost 
over time until we reach the long-time residual, Ii. The 
contributions that connect different states at the same 
energy are lost when we reduce down to the Mazur con- 
tribution directly described by the conservation laws, l2, 
but all the contributions from different c's still contribute. 
We find only the average contribution over the degener- 
acy label. The final reduction to the Hamiltonian-only 
Mazur contribution involves an average over all the ma- 
trix elements at each energy, I3. 

We can now answer the question of when the commut- 
ing algebra of conservation laws is enough to describe 
the long-time residual of the current: Formally, the long- 
time residual is exhausted by the conservation laws only 



when the current operator restricted to each energy sub- 
space is the identity in the basis which diagonalises the 
symmetries. Given that there is a conservation law which 
exhausts the long-time residual, jc, we still have the issue 
of when we manage to capture it in an arbitrary collec- 
tion of conservation laws. In the original argument Q a 
"maximal set of mutually commuting conservation laws" 
is used as a starting point. This is not suffieient! In 
general the 'mutually commuting' is a technical problem: 
Note that saturation fails when the relevant conservation 
laws are non-Abelian, such as 



2) (2 



1)(2 



2)(1 



(38) 



where only the first two are maximally mutually commut- 
ing, but the current operator could include some of the 
1)(2 I + I 2)(1 I conservation law, which is inaccessible. 
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The degeneracy is clearly a crucial issue, since if the 
eigenbasis is non-degenerate then the labels c and m are 
irrelevant and we find that 



^ ^ jnnjnn i 



(39) 



and the Hamiltonian itself generates the long-time resid- 
ual and the conserved piece of the current is an elemen- 
tary function of the Hamiltonian 



Jc 



I n)jnn{n \=^jnnP7i[H], 



with 



Pn[H] - n 



H 



(40) 



(41) 



the relevant projection operator. Indeed, for this case the 
Hamiltonian generates the only conservation laws and the 
system should not be thought of as being integrable in 
any way similar to the classical analogues at all. 

For any mutually commuting set of conservation laws 
we can construct a single conservation law from which all 
others may be constructed. We may choose 



Ac = ^^\ n,c, ■m)Anc{n, c, m |= ^ AncPnc, (42) 



nc m 



as in equation ()28|1 . and if all the Anc are distinct then 

{Ac — A„'c') 



and consequently 

Ca 



nc 



(43) 



(44) 



is an elementary function of Ac- In a calculation of a 
Mazur inequality, in principle, we can replace a set of cho- 
sen conservation laws with powers of a single law which is 
an arbitrary linear superposition of the original, because 
such a law is generically non-degenerate with respect to 
the chosen set. 

In classical problems we think in terms of phase-space 
and trajectories: Each conservation law corresponds to a 
restriction on the dimensionality of the subspace that a 
trajectory is moving in, with each law corresponding to 
the loss of a dimension and a complete set of N conser- 
vation laws in 2A'^-diniensional phase-space leading to an 
A^-torus remaining with the trajectories cycling endlessly 
about the N handles of the torus at independent rates, 
for an integrable system. In quantum problems the ana- 
logue is that an eigenvalue of an operator is fixed and the 
system is restricted to the projected subspace for which 
the eigenvalue is fixed. The degeneracy in that subspace 
amounts to the residual freedom that remains to the sys- 
tem. Any additional conservation laws which commute 



with the Hamiltonian can be simultaneously diagonalised 
leading to a reduction in the projected subspace that the 
system can be restricted to. This analogy very much 
follows our development and we have analysed our long- 
time residuals in terms of these restrictions. The natural 
quantum analogue to a classical integrable system is then 
one for which the simultaneous diagonalisation of all the 
conservation laws leads to a unique basis and no residual 
freedom. 

For an integrable system we anticipate that the la- 
bels n and c are sufficient and that the label m becomes 
redundant. It is still not clear whether the conserva- 
tion laws exhaust the long-time residual, because the 
current operator could connect states at the same en- 
ergy but with different values of the other conservation 
laws. This eventuality would correspond directly to the 
previous comment that even a complete set of mutually 
commuting laws docs not guarantee to exhaust all non- 
Abelian conservation laws. We now move on to a particu- 
lar case, and analyse these formal ideas for a concrete ex- 
ample in an attempt to resolve the failure of conservation 
laws topredict the long-time residual of the Heisenberg 
modelHfnl. 



III. LONG-TIME RESIDUALS AND 
CONSERVATION LAWS IN THE XXZ MODEL 

The particular system that we use for the current nu- 
merical investigation is the XXZ model 



_ (45) 
which is integrable [iJI- For this model the conservation 
laws have been established^. The appropriate boost 
operator is 



(46) 



and a sequence of conservation laws can be generated by 



Cn+l 



B, Cn 



(47) 



with Co = H, in notation consistent with our formal 
analysis. These conservation laws are local, in the sense 
that they arc a sum of translated copies of interactions 
which involve a finite contiguous set of sites. Indeed, each 
subsequent law involves an extension of the length of the 
interacting region by a single site. For all the undeniable 
mathematical beauty of this construction, unfortunately, 
it is of no practical use to us! To investigate saturation, 
we have been forced into analysing finite systems where 
the number of conservation laws is controllable and the 
extensive nature of the boost operator makes it incom- 
patible with the periodic boundary conditions necessary 
to our calculations. 
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We have employed a more cumbersome approach to where 
creating the conservation laws. The construction involves 
the 'R-matrix' 



A = cos (5 < 1 
A 1 
A = cosh S > I 



S{X) = 



sin 1(1 + A) 



sm ■ 



S{X) 



5(A) = 1 + A 
sinh|(l + A) 
sinh I 



C(A) = 



cos 1(1 



A) 



cos ■ 



C(A) 



C(A) = 1 

cosh|(l + A) 
cosh I 



(49) 



in terms of which one can construct the monodromy 
matrixQ for N sites with periodic boundary conditions 



f[A] =tro i?oi[A]i?o2[A]...i?ow[A] 



(50) 



The trace is over the spin So which is a dummy spin. This 
matrix controls the conservation laws using the result 



f[A],f[A'] 



0, 



(51) 



for all A, A' gR, which guarantees the linear indepen- 
dence of any operators that can be constructed from this 
matrix. The local conservation laws are derived from 



In (f[X]f[0]-^ 



(52) 



with Co = H, and all subsequent operators being inde- 
pendent. 

The previous formal ideas are straightforward to anal- 
yse on a finite system: Employing a 'dummy' spin to 
describe So it is elementary to calculate the monodromy 
matrix as a function of both A and i5. This real but non- 
Hcrmitian matrix may be diagonalised and is almost in- 
variably (respecting ~ J2i ^f) non-degenerate. This 
provides an essentially unique basis in which all the con- 
servation laws are diagonal. The final degeneracy label in 
our analysis becomes irrelevant and our analysis reduces 
to 



In = i Vp-/3^" V f',/," 

-^^> 2" / ^ / V Jnn' Jn' ni 

n n'cc' 

-^i ^ / J / J ■JnnJnn i 

n cc' 

/, = - V e--^^" V i'^" 

2" / ^ / J ■JnnJnn') 

n c 

h = -Y^ e-'^'" — V 7^^ V 7^'^' 

^ / ^ ^ / ^ Jnn / ^ Jnn i 



(53) 
(54) 
(55) 
(56) 



Our formal considerations involved an arbitrary oper- 
ator, j, which we designated a 'current'. For the XXZ- 
model there is a physical aspect to the current, which can 
be constructed from either an infinitesimal 'twist' in the 
boundary conditions (corresponding to a magnetic field 
through the loop) or the local continuity equation 



— (x)+V.r(x)=0 



1 

ih 



H, 



1 



Ax 



J. 



= 0, 



which is solved by 



J 



(57) 



(58) 



Up to a constant, this provides our chosen current oper- 
ator 



_ r, / dy cx QX dy 



3 



(59) 



We note that the Hamiltonian contribution vanishes for 
this operator. This can be proved using the symmetry 



U = exp 



(60) 



a rotation through tt about the x-axis, which preserves 
the Hamiltonian but reverses the sign of j. In the basis 
that simultaneously diagonalises H and U we find that j 
is necessarily off-diagonal and consequently I3 vanishes. 
Note that U does not respect but does respect the 
symmetries C° and so does not impinge on the standard 
conservation laws. Although in principle Ii and I2 are 
distinct, for this model we find that the two quantities 
agree exactly and in the unique natural basis the relevant 
matrix elements of the current arc diagonal. 
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As yet we have not encountered our conundrum: If 
we employ all the conservation laws then in our formal 
sense we achieve all the long-time residual for the cur- 
rent model. It is not straightforward, in practice, to em- 
ploy the eigenbasis and one would like to employ a single 
conservation law that has a non-trivial overlap with the 
current to demonstrate whether the long-time residual 
exists at all, let alone is being exhausted. It is this issue 
that has proven problematic for the current model, where 
at 'half-fiUing' (corresponding to the subspace S^=0) it 
has been shown that none of the conservation laws can 
contribute to the long-time residual, essentially because 
they respect U for this filling and then any expectation 
value involving both the current and any combination 
of the local conservation laws vanishes, eliminating the 
contribution to the Mazur inequality. 

So far we have discussed the pure XXZ model, but 
there is a natural extension to include a magnetic field 

H = Ha-BY,S:, (61) 

i 

which activates the local conservation laws and allows a 
non-trivial contribution to the long-time residual to be 
obtained from the simplest lawfl^. This additional term 
commutes with the original Hamiltonian and all the local 
conservation laws and completes the uniqueness of the 
diagonalisation. Note that U is no longer compatible and 
we have to use spatial inversion symmetry to demonstrate 
that the Mazur contribution from the Hamiltonian now 
vanishes. We employ this field in our calculations as a 
natural parameter to contrast the different behaviours 
and assess whether or not half-filling is a pathological 
case. 

We are now led to a crucial issue, that of the compati- 
bility of the current with the conservation laws: We have 
chosen to complete the intcgrability of our model using 
the z-component of total-spin (which would be forced in 
the presence of the field). If, instead, we chose to com- 
plete the isotropic case, A=l, with the x-component of 
the total-spin (or used U itself for the general case) then 
U remains a symmetry and the current operator is neces- 
sarily off-diagonal. For this case there is no Mazur contri- 
bution and the current operator is incompatible with the 
maximal set of mutually commuting conservation laws. 
Recognising that this non-Abelian aspect to the conser- 
vation laws can be critical, we now return to the case 
where we employ the z-component of total-spin and the 
current is compatible with the conservation laws. 

The monodromy matrix provides us with a conserva- 
tion law for each value of A, and by varying A we have 
access to a set of independent conservation laws. Indeed, 
for a chain of length N we have N linearly independent 
conservation laws. Unfortunately, these are not the pure 
local conservation laws but for each order an additional 
polynomial of lower order local conservation laws is in- 
corporated into these conservation laws. In the issue of 
saturation this does not cause problems, but making par- 
ticular deductions about the local conservation laws is 



usually made impractical. Although we generate a 'com- 
plete' set of conservation laws there is an important sub- 
tlety. A particular conservation law can be used to gener- 
ate more: Any polynomial combination of a conservation 
law is another such law and as we have seen, the pro- 
jection operators, which project onto eigenspaces, are a 
complete set of conservation laws which can be generated 
by any particular conservation law. Further, it should be 
appreciated that for N sites we have a state-space which 
scales as 2^ and is exponential in system size. Usually, 
the number of derived conservation laws is exponential 
in the system size and the total number of conservation 
laws required to get a complete set is disturbingly large 
even for quite a small system. 

We perform a sequence of numerical calculations: For 
a chain of length TV we choose N distinct values of A and 
create N distinct monodromy matrices. We then calcu- 
late our four quantities as a function of applied field for 
the set of all N such conservation laws. Note that we 
choose to orthogonalise these conservation laws because 
of the additional stability this offers for the case when the 
conservation laws arc over- complete. We then create new 
conservation laws by taking the N{N+\)/2 distinct prod- 
ucts of these monodromy matrices and include these in a 
larger calculation with both linear and quadratic combi- 
nations and calculate a second estimate to the Mazur in- 
equality. We then calculate yet further conservation laws 
by taking the N{N + 1){N + 2)/6 distinct triple prod- 
ucts of these monodromy matrices and including them 
into a third Mazur estimate. We took this procedure 
to quartic order but proceeded no further. The results 
for a couple of representative cases, one metallic and one 
insulating, are provided in figure 1. The conundrum is 
clearly observable: The contribution to the Mazur in- 
equality from the conservation laws smoothly vanishes at 
zero field and some of the long-time residual is lost. For 
a finite system we can exhaust the attainable conserva- 
tion laws and this observed loss is genuine! Fortunately, 
the resolution of this issue is straightforward: One im- 
portant conservation law was missing from our analysis, 

, the z-component of total-spin. The monodromy ma- 
trix is only non-degenerate when wc additionally respect 
this symmetry too. Including this conservation law into 
our previous calculations, by including times the pre- 
vious conservation laws as 'independent' laws, provides 
figure 2. The problem at 'half-filling' has been alleviated 
and for small systems all the conservation laws can be 
exhausted and we can achieve the full long-time residual. 

For the particular case of the XXZ model and the stan- 
dard spin current, we now have an answer to our initial 
question: an initial current does decay as a function of 
time, except at A=0 or oo, and the long-time residual 
is exhausted by the known conservation laws. In order 
to exhaust this residual one must include into the 
list of conservation laws, because without this law the 
residual is not exhausted and indeed at half-filling the 
local conservation laws offer nothing. Our results are for 
finite systems and there arc definite problems in the thcr- 
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FIG. 1: Long-time residuals and estimates to the Mazur in- 
equality. We offer the original current, the exact long-time 
residual and a sequence of three estimates coming from the 
conservation laws. The third calculation exhausts all the 
available conservation laws. N=8, 5—1, /3=1 and (a) A <1 
(b) A >1. 

modynamic limit. It is 'known' [lol [l]|. from the Bethe 
Ansatz, that at half-filling the long-time residual shows 
interesting behaviour: For A <1 we have a finite resid- 
ual that smoothly vanishes as A (—>1 and vanishes for 
A >1. A single local conservation law must find diffi- 
culty in providing this behaviour and it is our need for 
the combination of 5^ with the local conservation laws 
which overcomes this difficulty. 

Since any particular value of A generically provides a 
non-degenerate monodromy matrix, we could have used 
a single such conservation law to exhaust the conserved 
current. We have demonstrated that this is so on very 
small systems, but the similarity of the laws provided by 
powers of this one law lead to unpleasant instability is- 
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FIG. 2: Long-time residuals and Mazur estimates. We offer 
the original current, the exact long-time residual and a se- 
quence of four estimates coming from the conservation laws, 
both with and without the total-spin 5*^ included. The fourth 
calculation exhausts all the available conservation laws. N=8, 
(5=1, 13=1 and (a) A <1 (b) A >1. 

sues numerically, which, together with the added control 
of completeness issues, made our chosen procedure more 
effective. 

Since we are dealing with the physical current in the 
system, but the fundamental theory should be applicable 
to an arbitrary operator, one can ask if the results are 
in any way special for the physical current. We repeated 
our calculations for the longer-range 'currents' 

jn = ( — «) ^ i^^t^j+l+n ^ + ; (62) 

j 

where jo is the physical current and ji and j2 are ac- 
cessible to our available system sizes and they exhibit 
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essentially no new physics at all: The operator tj is still 
available to show that the local conservation laws offer 
nothing in the absence of a field and the analogue plots 
to figure 2 are equivalent in all respects. 

What can we say in general? At the formal level 
our results are clear: We see four natural quantities, 
^0 > ^1 > -^2 > -^3, and in principle they are all dis- 
tinct. The initial current, controlled by /q, only remains 
undiminished if the current operator is a conservation law 
itself, as happens for the case A=0. The long-time resid- 
ual, controlled by /i, may be considered to come from 
conservation laws, but only if we allow the inclusion of 
non-Abelian laws. If we have a set of mutually commut- 
ing conservation laws, then we only exhaust the long- 
time residual when the previously considered conserved 
current is the identity in the basis which diagonalises the 
conservation laws and also conserves any residual degen- 
eracy. 

We have shown that the symmetry at zero field causes 
the local conservation laws to become independent from 
the conserved current and we have resolved this by in- 
cluding the z-component of total-spin to successfully rep- 
resent the entire conserved current, but all this was for 
finite systems: What happens in the thermodynamic 
limit? To exhaust the current conservation law we need 
to generate the long-time residual to machine accuracy 
from a Mazur calculation and this we can achieve with 
all accessible systems which do not have to be large to 
be convincing. To analyse the thermodynamic limit we 
need to finite-size scale and our relatively small systems 
are only indicative and not at all convincing. Never- 
theless, we believe that in the thermodynamic limit the 
point at zero field becomes singular and that at all other 
fields the local conservation laws exhaust the conserved 
current. One analytical consideration is that in the ther- 
modynamic limit we should get the same answer if we 
canonically fix the magnetisation or grand-canonically 
fix the field: Since when the magnetisation is fixed there 
is no difference between the conservation laws with and 
without magnetisation (except perhaps when the mag- 
netisation vanishes), then there can only be one quantity 
generated by these laws, the full conserved current. We 
have incorporated the canonical calculations into figure 
2 and it is plausible that they converge to the long-time 
residual. We have also finite-size scaled the fraction of 
the long-time residual provided by the conservation laws 
in figure 3 and it is easiest to believe that in the ther- 
modynamic limit the point at zero magnetisation is an 
isolated puncture. 

Although we appear to have a large number of degrees 
of freedom in our calculations (4096 for iV=12 for ex- 
ample) a brief analysis of the available degeneracy might 
suggest that our systems are too small to be relevant: 
Once we have extracted the 5*^ degeneracy employing 
a field, although at first sight there appears to be a 
huge residual degeneracy, in fact this degeneracy is al- 
most all lifted by inversion symmetry and the further 
residual degeneracy is woefully small! Indeed, we may 



Convergence 




Field 



FIG. 3: Overlay of the ratios of exhausted Mazur estimates 
without the total-spin included. A'^=4,6,8, (5=1, /3=1 and 
A <1. 

employ the first non-trivial local conservation law, the 
energy current [T^ . to lift the inversion symmetry and up 
to A^=ll this fully diagonalises the system. Additional 
degeneracy emerges for iV=12, but this is much less than 
is naively expected: We might expect that at order N 
the N distinct conservation laws were functionally in- 
dependent, but for systems up to 7V=11 all N laws are 
describable as functions of the Hamiltonian, the energy 
current and the z-componcnt of total-spin! The role of 
degeneracy in the thermodynamic limit is currently not 
accessible to computers, even at the level of a hint! 

For the Heisenberg model, we have established that for 
a finite system we can expect that 

(j(^).?(0)) = (jJc), (63) 

with 

j,=T[S\C''] (64) 

for some function of the relevant conserved quantities. 
The final issue is whether we can use this knowledge to 
prove the existence of a non-trivial long-time residual in 
the thermodynamic limit. Employing the first non-trivial 
conservation law, the heat current, it was shown that 
there is a non-vanishing residual at finite temperature 
for all fillings except half-filling p^. As we have seen, at 
half-filling we also require to use the conservation law 
5^, so it seems plausible that we might be able to use a 
combination of both the heat current and the total spin 
to find the residual at half-filling. Unfortunately, this 
does not appear to be true and a verification of the long- 
time residual at half-filling using the Mazur inequality 
must await further investigation! 

A natural technique is to Taylor expand the general 
conservation law and create a trivial conservation law 
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from which the long-time residual can be deduced from 
the Mazur inequality 



(j(oo)j(O)) > 



(C'tj(0))(j(0)g) 
(CtC) 



(65) 



Clearly all the (7° have the wrong symmetry and itself 
also has the wrong symmetry and so we might like to start 
from the law 



C = S'^Ci, 



(66) 



for example. We can analyse this proposal in detail for 
the limit A=0 where the Jordan- Wigner transformation 
reveals the solution as a non-interacting free-electron gas. 
Employing 



3 = 2^sinfc/^^/fe, 



Ci^^sin2fc/t/fe, 



(67) 



flh 



we can verify that 

(jC^S^) = 2^sinfcsin2fcnfc[l-nfc][l-2nfe], (68) 



with 



nk 



1 + eP^" 



(69) 



at half-filling and does not vanish, as expected. Unfor- 
tunately, in the thermodynamic limit this conservation 
law offers nothing because it scales with the length of 
the chain but 



l\2i QZ\2\ 



z\2\ 



(70) 



scales with the square of the chain length rendering the 
contribution negligible! Note that the case A=:0 is spe- 
cial: There is a second ladder of conservation lawsQ 
which includes the current operator itself. When A ^0 
the conservation laws 



C2m = ^ cos(2m + l)kflfk, 

k 

C2m-i = ^ sm2mkflfk 



(71) 



smoothly connect to the local conservation laws whilst 
the laws 



D27n = ^cos2mfc/2/fe, 

A; 

-D2m-i = ^ sin(2m - l)kflfk 



(72) 



are broken (except Dq which is not part of the ladder). 
Noting the symmetries of these classes, the law 



C = S^ 



CgCq 



(73) 



contributes from the appropriate sector and yields a fi- 
nite Mazur contribution. This law arises from a linear 
superposition of the local conservation laws but it itself 
is clearly long-range. 

One of the central current issues is that of finding a 
conservation law for which a single application of the 
Mazur inequality would provide incontrovertible proof of 
the existence of a long-time residual. We already know 
that this must be non-trivial for the current model, since 
the long-time residual is known to vanish for A >1 and 
hence any prospective conservation law must respect this 
fact, a very non-analytic prospect. If we consider the gen- 
eration of part of the long-time residual using a single 
conservation law, then we encounter an immediate prob- 
lem if we attempt to employ local laws. If we consider 
the conservation law 



S'C, 



(74) 



where we require to use a local conservation law and the 
z-component of total-spin, then we may replace this with 



S'-{S') C - {C) -{ - {S') C- (C) ), (75) 



since neither operator has an overlap with the current on 
its own. Now the Mazur inequality tells us that 



(.Koo)j(o)) ^ ^^'0 (^^ ' ^^^0 ' ' ' 

6-(o)j(o)) - (^--(5^)) [c-{c)) {c-{c)y 



(76) 
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and note that, employing notation 

6A = A-{A), 

we have 



Finite-size Scaling 



(77) 



{5S'5S^5C6C) = {5S^^5S~-){5C5C) 

+{\5S''5S' - {SS^'SS')] UCSC ~ {SCSC)]), (78) 

where the first term is likely to dominate. These ideas 
are essentially true whatever the operator C, but if C is 
local and the correlations in the system are short-range 
then we can deduce something from this: If we have that 



A 



for all such local operators then 



(79) 



(80) 



and if i and j become far apart we might expect the cor- 
relations to become irrelevant. This correlation function 
might then be expected to scale linearly with system size. 
Similarly, including a C then 



{6A5B5C) ^ ^{5A,5Bj5Ck) 



(81) 



ijk 



also might be expected to scale linearly with system size. 
The consequence of this is that any local contribution to 
the Mazur inequality would be expected to vanish in the 
thermodynamic limit. In practice, the correlations are 
likely to be power-law in nature, but the physical idea 
behind the failure remains. 

The hunt for an elementary conservation law that will 
provide a non- vanishing Mazur contribution at zero field 
probably necessitates non-local conservation laws, but we 
should remember that the monodromy matrix expanded 
in the manner that we employ does yield non-local lawsQ 
and further we are complete for our finite systems so we 
employ an arbitrary linear superposition. Based on the 
previous limit of A i-^O we might suspect that it is only 
an arbitrary superposition over all lengths that is cru- 
cial and not the non-linear combinations of laws. Our 
final calculation offers a test of this idea through a finite- 
size scaling calculation of the fraction of the long-time 
residual provided by the linear conservation laws, the 
quadratic conservation laws and so on, order by order, in 
figure 4, and unfortunately the contribution scales away, 
reinforcing the fear that a truly non-analytic function of 
the conservation laws is required. 



IV. CONCLUSIONS 

Integrable systems can have anomalous conductiv- 
ity, with currents flowing indefinitely once started. As 




FIG. 4: Finite-size scaling plot of the fraction of the long-time 
residual achieved by the linear contribution to the conserva- 
tion laws for the case A'^=4,6,8,10, 5=1, P=l and A <1. 



suggested by Suzuki Q, this long-time residual can be 
thought of as some part of the current operator actually 
being one of the conservation laws of the system. We 
agree with this statement but refine one of the issues 
arising: Given a set of conservation laws, when is the 
conserved current either partially or wholly described by 
them. One might anticipate that a 'complete set of mu- 
tually commuting' conservation laws would generate all 
possible conservation laws, but this is incorrect: In gen- 
eral, additional independent conservation laws can ex- 
ist with non-trivial commutation with the complete set. 
This is noi just a formal issue, as can be seen by example: 
The isotropic Heisenberg model with the standard spin 
current. The local conservation laws do not completely 
diagonalise the system and there are a variety of ways of 
completing the conservation laws. We can simply diag- 
onalise one of the components of total-spin as well, for 
example. If we choose to use the z-component of total- 
spin then we can generate the full conserved current but 
if we choose either the x- or y-component of total spin 
then we can generate no long-time residual whatsoever! 
The current must be compatible with the complete set of 
mutually commuting conservation laws. 

Using numerics on finite systems, we show how to gen- 
erate the current conservation law for the XXZ model: 
If we combine the local conservation lawsQ with the z- 
component of total spin then we can exhaust the long- 
time residual and consequently form the complete con- 
served current. Since the local conservation laws yield 
no long-time residual at all, in the absence of a field, any 
conservation law which has an overlap with the current 
must involve both the local conservation laws and the 
z-component of total-spin. 

The final issue addressed is as to the reason for the fail- 
ure of any elementary conservation laws to provide even 
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part of the conserved current for the case with no field. 
We observe that since we require to employ both a local 
conservation law and the ^-component of total-spin then 
for short-range interactions we need a coincidence of three 
local operators to yield a contribution and a consequent 
loss of this contribution in the thermodynamic limit. It 
is not easy to use the conservation laws to demonstrate 
a long-time residual in practice! 
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